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CARDIOVASCULAR MODEL^^ 
(lumped) 



Multiscale model for VHP research 





CEREBROVASCULAR MODEL ^ 
(lumped) 


EYE MODEL (lumped) 


EYE MODEL (finite eiement) 



http://www. alpha 
m/tcd-cai 
doppler 




http://grdedev.ferris.edu/~oddenj/GRDE%20228/how-it- 

works/Diagrams.html 


- Linningeretal. (2009) 


CNS model includes intra/extracranial 
cerebrospinal fluid (CSF) and cranial blood 
compartments 

For details on other modules, see companion 
works for IWS2015 by Ethier et al., Feola et 
al., Nelson et al., and Price et al. 
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CNS Blood flow and pressure model 





Several lumped CNS models exist. Our 
starting point was a model that had been 
applied to microgravity (pg) (Stevens et al., 
2005; Lakin et al, 2007); 

Time-dependent model composed of 6 fluid 
compartments (nodes) 

• 3 vascular: 

• Intracranial Arteries (1 ) 

• Capillaries (2) 

• Venous Sinous (3) 

• 2 cerebrospinal fluid 

• Ventricular CSF (4) 

• Extraventricular CSF (6) 

• 1 Brain node (5) 

Boundary conditions at cranium and whole- 
body interaction provided by extracranial 
nodes 

• Central Arteries [A] 

• Central Veins [V] 

• Thoracic Space [Y] 


Q = Flowrates between compartments (ml/min) 
C = Compartment compliance 


Brain (5) | 

^ 1 Qo, 


Intercranial 

Arteries 


Central 

Arteries 


Capillaries (2) 


Ventricular 


CSF (4) Q, 



Brain (5) 


\ Venous^ 

^ 63 ) ( C53 


Extraventricular CSF (6) '^6V ) Central 




Veins 


Thoracic Space [Y] 


- Stevens et al. (2005) 
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Governing Equations 




• Defining the pressures in the 6 compartments 
as dependent variables, the system is modeled 
in matrix form as a system of ordinary 
differential equations: 


C [dF^ + ZP = S 
_dt 



■^Ai{Pa ■ ^A|Sin0) - Qi2 


■ Qi 4 ■ 1 ^ 23^25 (^2 " 


Zjy{Py - G^i^sin9) 


■ 1 ^ 25^25 [^2 ■ 

C6y{dpy/dt) + 
Qv(dpv/dt} + Zftv(Pv) 


C compliance 
G gravity 

K filtration coefficient 
P pressure 
Q flow rate 

S source/forcing terms 


Z fluidity ~ 1 /resistance 
0 tilt angle 
n osmotic pressure 
a reflection coefficient 


Note that G is explicitly included in the forcing 
terms in S 


National Aeronautics and Space Administration 


Glenn Research Center 4 









MATLAB Implementation 
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Verification Tests 



* 20 independent verification tests that included variation 
in hydrostatic pressure 

• 3 independent users of the code 


TEST 

Short-term head down tilt (HOT) 
Long-term HOT 
Microgravity 

Blood-brain barrier influence 

Verification tests also had a validation component 

• Used Lakin and Stevens equation structure and parameters, but 

• Developed independent implementation, arterial pressure that drives 
unsteady response and solution methodology 
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Short-term head down tilt 


• Tests called for monitoring of changes in pressure differences pre- to post-tilt: 



Tilt 

angle 

n 


-6 


Expt 

[ 1 ] 


A(Ps-Pv) (mmHg) 


Model 

[2] 


-10 


3.1 


3.86 


-15 


Our 

Expt 

Model 

[3] 


3.3 

3.66 



A ICP(mmHg) 


Model Our 

[2] Model 


2.10 to 3.70 2.20 


6.1 


5.18 to 7.78 


5.46 


Conclusion: The current model agrees with prior 
experimental and numerical work. 


[1] Katkov and Chestukhin 
(1980) 

[2] Stevens et al. (2005) 

[3] Murthy et al. (1992) 
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Long-term HDT and Microgravity 







Time (s) 

25. Capillaries to Brain 



Time (s) 

x^ txtraVentricular CSF to Veins 



Time (s) 



o 

S 


7J 

>. 


m 

</> 


Condition 

A iCP (mmHg) 

Modei [1] 

Our Modei 

Long-term HDT 

4.9 

4.9 

rig 

<0 

-3.1 


[ 1] Stevens 
et al. (2005) 


Conclusion: Using their parameters, our predicted AlCP is 
consistent with the prior model in jug and long-duration HDT, 
. but are their parameters correct? 
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Blood-brain barrier influence 


• Later work by the Stevens/Lakin team hypothesized that the blood/brain 
barrier might weaken in |Lig 

• In Lakin et al. (2007), they performed a sensitivity study for a hypothetical 
change in the reflection and filtration coefficients 

• This changed their findings on ICP in jug 


ICP isobars, Lakin et al. (2007) 



13 mmHg 

15 mmHg 

17 mmHg 

19 mmHg 

O TestO 

• Testl 
O Test2 
O Test3 

• Test4 


TestO 

Testl 

Test2 

Tests 

Test4 

sigma 

K 

Simulated 
ICP (mmHg) 

Target 
ICP (mmHg) 

Error (%) 

1.000 

0.066 

13.78 



0.583 

0.052 

15.15 

15 

1 

0.665 

0.105 

17.18 

17 

1.06 

1.081 

0.064 

13.23 

13 

1.77 

0.438 

0.113 

19.14 

19 

0.74 


Conciusion: Our model agrees 
with literature results to within 2% 
or better. 


Revising prior findings, authors concluded that ICP could increase in |ig. 

But how do we assess the credibility of this claim? 
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Preparing for |iig simulations 


• Before weighing in on the potential change in ICP in \ig, we need to: 

• Re-assess parameters used by Lakin/Stevens based on the most current 
VHP research 

• Quantify uncertainty in model parameters 

• Define a physiological envelope for parameters that will be relevant for 
the astronaut corps on orbit 

• Perform sensitivity studies over a much larger parameter space 

• Examine model predictions against independent studies in HDT, |Lig, and 
postural change, particularly for chronic conditions. We need our model 
to do a good job in predicting: 

• Volumes of intra/extracranial CSF compartments 

• Volumes of intracranial blood compartments 

• Only after these steps are taken can we make intelligent predictions 
about |LXg response 
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Sensitivity analysis 


• We are analyzing this system by testing model sensitivity 

• Parameters include: compliances, resistances and filtration coefficients 

• Each described by statistical parameters 

• Mean and range of variation (variance) 

• Distribution of variation (density function) 

• Methodology 

• Partial Rank Correlation Coefficient (PRCC) Analysis 

• Provides the linear relationships between two variables 

• one input parameter and one output parameter 

• All linear effects of other variables are removed after rank 
transformation 

• Rank Transformation: transforms nonlinear monotonic relations to linear 

• Latin Hypercube sampling 

• Efficient method to randomly characterize the sets of combined parameters 

• Many independent runs with randomly chosen parameter sets provide 
statistics on the system response 
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Conclusions 




• A CNS lumped parameter model has been produced based on the 
model developed by Lakin and Stevens 

- Our solution methodology and computational platform is unique 

• Our model has been tested and verified 

- ICP predictions agree with Lakin/Stevens in 20 cases of acute and chronic 
|Lig and HDT 

• CNS model infrastructure is complete, but additional work is needed 

- Re-assess parameters used by Lakin/Stevens 

- Define flight and flight analog derived parameter ranges 

- Perform parameter sensitivity studies 

- Validate against the latest VHP research 

• In the future this model will be 

- integrated with lumped CVS and eye models 

- Used to establish spaceflight responses with fidelity sufficient to supply 
boundary conditions for more complex VHP eye simulations. 
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